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^ | ■ For Lindblad's master equation of open quantum systems with a general 

■ quadratic form of the Hamiltonian, the propagator of the density matrix is 

analytically calculated by using path integral techniques. The time-dependent 
density matrix is applied to nuclear barrier penetration in heavy ion collisions 
I/"") ' with inverted oscillator and double- well potentials. The quantum mechanical 

decoherence of pairs of phase space histories in the propagator is studied and 
shown that the decoherence depends crucially on the transport coefficients. 
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I. INTRODUCTION 



In many problems of nuclear physics and quantum optics, where one deals with open 
quantum systems, the memory time of the environment is very short and a Markovian 
approximation is suitable. Disregarding the averaging over the intrinsic degrees of freedom, 
one can consider the open system starting from the general Markovian master equation for 
the reduced density matrix of the collective degrees of freedom as given by Lindblad |l 



Ho, p] + ^-E (ftA Vj + ] + pVj + }) ■ (!) 



dp(t) i 

dt = ~n v "^ ri ' 2h 

Here, Hq is the Hamiltonian of the collective subsystem and Vj are operators acting in the 
Hilbert space of the subsystem. The terms in the sum of Eq. (JiJ) are responsible for the 
friction and diffusion and supply the irreversibility of the dynamics of the open quantum 
system. Omitting these terms we get a standard form for the evolution equation of the 
density matrix for closed systems. This equation and similar equations were used e.g. in 
Refs. HH. 



Path integral methods are a conventional tool to describe open quantum systems flO|-|T8 
Here, we use results of Strunz 0] elaborated with path integral techniques and derive 
analytical expressions for the time- dependent density matrix for a general Hamiltonian of 
quadratic form with an inverse oscillator potential which can be applied to the description 
of fission and fusion through potential barriers in nuclear physics. The decoherence of pairs 
of phase space trajectories will be studied in the semiclassical limit for different choices of 



the effects of the environment on the system. As was shown in JTTJ, the initial Gaussian 



distribution remains to be Gaussian in an oscillator potential. We extend this statement 
for any quadratic form of the Hamiltonian of the subsystem. By a direct numerical solution 
of Eq.([IJ) we consider the evolution of the density matrix in time in a double- well potential 
under various sets of transport coefficients. Such potentials are more useful and realistic to 
investigate nuclear fission problems than inverse oscillator potentials. 

II. PATH INTEGRAL PROPAGATOR AND DECOHERENCE 

With the propagator G(q,q',t;q Ol q' o ,0) of the density matrix one can find the density 
matrix < q\p(t)\q' > (in coordinate representation) at any time from the initial one < q\p(t = 
0)\q>>: 

<q\p(t)\q'>= J dq J dq' G(q,q',t;q ,q' ,0) <q \p(t = 0)\q' > . (2) 

In the one-dimensional case an expression for the phase space path integral of the propagator 
corresponding to ([I]) was derived in [jlO 



as 



(ff,<) (</',*) 



G(q,q',t;q ,q' ,0) = J D[a) J D[a'\ exp iiS[oi\ a')), 

(90,0) (^,0) 
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S[a; a'} = S[q,p; q',p'} = J rfr{g(r)p(r) - # e// (g(r),p(r))} 

o 

dT{q'(T)p'(T)-H: ff ( q '(r),p\r))} 



o 

t 



- iW dT{VMr),p(r))V;(q'(r),p'(T))}, (3) 

3 

with phase space paths [a] = [q, p] , where q and p are the position and momentum, respec- 
tively. The effective Hamiltonian is given by 

H eff =H - \Vj\ 2 . 

Z 3 

Here, the quantities H , \Vj\ 2 , Vj and V* are the Wigner transforms of the operators H , 
V^Vj, Vj and in (1), respectively. 

Choosing an inverse oscillator potential, we write the Hamiltonian of the collective sub- 
system in a more general quadratic form 

= 2m ~ ~lT q + f (W + qp) ■ (4) 

The environment operators are assumed as linear 

V j = A j p + B j q, V+ = A*p + B*q, i = 1,2. (5) 

As shown in [|lOj| for the similar case of an harmonic oscillator, the integrals in (§) over the 
momentum yield Gaussian integrals and can be evaluated. Then the propagator is reduced 
to path integrals in coordinate space [[l(J : 

9(0 <?'(*) 

G(q,q',t;q ,q' ,0)= J D[q] J D[q') exp {^S[q; <f\), 



S[q\ q'} = S d [q] - S d [q'} - iHXt + <%, q'} + ~D[q, qf. (6) 

In Eq. (|6|) the classical action of the isolated system S c i, the phase function <&[q, q'] and the 
square of the decohering amplitude D[q, q'\ can be expressed as 

jdr^me + ^V], (7) 







* A - A r A A * 



dr(qq' — qq') + m-^—^ — - J dT(q'q' — qq) — m p 9 J dr(q 2 — q' ), (8) 
o oo 

D[q, qf = I | (D pp + m 2 \ 2 p D qg + 2m\ p D pq ) J dr(q - q') 2 



2m(D pq + mXpDgg) J dr(q - q) (q - q) + m 2 D qq J dr(q - q') 2 I . (9) 
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The quantum mechanical diffusion coefficients are D 



<ii 



Y^\AA 2 for coordinate, D 



pp 



2^|-B.j| 2 f° r momentum and D qp = — -^KeJ2A*Bj for the mixed case. The frictional 

j j 
damping rate A = —lmJ2 A*Bj and the diffusion coefficients must satisfy the constraint: 



D D — D 2 



> A h /4 with D qq > 0, D pp > which secures the non- negativity of the 
density matrix at any time. The values A p = A + fi and X q = A — ji (A p + \ q = 2A) are fric- 
tion coefficients for coordinate and momentum, respectively. Both position and momentum 
undergo a direct damping and diffusion process in contrast to the classical case. If D[q,q'} 
increases with time for q ^ q', then the propagator suppresses the non-diagonal components 
of the density matrix. Thus, the interference between different positions q and q' becomes 
weaker. 

Since H e ff depends at most quadraticly on p and q, the path integrals are Gaussian. 
In that case a semiclassical solution of the path integrals with the method of stationary 
phases leads to an exact analytical evaluation of the propagator. First, equations of motion 
along the path trajectories q{r) and q'(r) (complex trajectories) are calculated with the 
condition of stationary phase 6S[q,p; q',p'] = with S of Eq. (3). The following equations 
for Qi = q + q', Q 2 = q — q', P\ = p + p' and P 2 = p — p' result, which are solved with the 
boundary conditions q = (q(0) = q , q(t), q'(0) = q ', q'{t)). 



fQi\ 
Pi 

\P2j 



f-X, 

muj 


V o 



m 








4iAn 
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A p 
mu) 2 
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-AiD P 
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X n 
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(Qi\ 
Pi 

\P2j 



(10) 



Next, the solutions q(r) and q'{r) of Eq . (|10D depending on the parameters q , q(t), q ' and 
q'{t) are inserted into the action function S[q; q'] of Eq.(6) and integrated over r. The square 
of the decohering amplitude is found as: 



D[q,q'Y = smh-^t}(A t (q 



%) 



A_ t (q(t) ~ q\t)f + B t (q - q' )(q(t) - q f (t))), (11) 



where ip — ycu 2 + (X q — X p ) 2 /A and 

Dpp - m(mu 2 D qq + (X q - Xp)D pq ) 



A t 



2h\ 



+ exp[2\t]ip 



2 Dp P - m(-2D pq X p + mD qq {uj 2 - 2AA p )) 



2h(u 2 — X p X q 
1 



2h\(X p X q - oo 2 ) 
sinh [2 ijjt] {Dpp + m(2D pq X p + mD qq (uo 2 - X p (X q - X p ))} 

cosh[2^]{-AD pp + m[(-2D pq + m(X q - X p )D qq )^ 2 



B t 



2h{uo 2 - Xp\ q ) 
X(-D pq (X q - A p ) + mD qq (u 2 + 0.5(A 9 - A p ) 2 ))]}, 

tl^-A \) {lp[ ~ Dpp + m ^ 2D ^ + mD ^ ~ 2AA ^ 
X[Dpp + m(2D pq X p + mD qq {u) 2 - A P (A 9 - A p )))]} 
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^sinh[(^ + X)t] 



{ip[-D pp + m{-2D pq \ + mD qq (u 2 - 2AA P ))] 



fi\{(jj 2 — XpX q ) 

+ \[D pp + m(2D pq \ p + mD qq (u 2 - X p (\ q - A p )))]}. 



(12) 



Similar analytical expressions are obtained for the action S c i and phase $[g, q']. Then, the 
propagator @ is finally evaluated as 



G(q, q\ t; q Q , q Q , 0) 



mip 



exp(At) ex P (iS R /h) exp(-Z% g']7(2£)) 



2irh sinh(^i) 
S R = S d [q]-S d [<A + $[q,(A 

{cosh^t - 4>)[ql - %}+ cosh(^t + <p)[q 2 - q ] 



2 sinh(^t) 

— 2cosh(0) cosh(\t)[qoq — q q'} — 2cosh(0) sinh(At) [q^q — q q]}, (13) 

where sinh0 = (X q — X p )/(2u). This propagator is correct for any quadratic Hamiltonian 
and is a generalization of the results of Refs. ||T0|-|12]1 where propagators were obtained for 
harmonic and inverted oscillators only For the initial density matrix (q(0) and p(0) are 
mean values) 



<q\p(0)W> = (2vra OT (0))- 1 / 2 

1 



x exp[ 



4M0) 



{(? - mr + - mn - =p{o)tf - <?)] 



n 1 



(14) 



the density matrix at time t is calculated with (0) and (13) as follows 



< q\p{t)W > 



: exp 



q + q' 



+ 



2na qq (t) 
icr pq {t) fq + q' 



2a qq (t) 



or in explicit form 

<q\p{t)\q> : 



?(*) (?-?') + (9-?') 



(9 " 

(15) 



2 sinh(V>t) 



exp (At) 



2tr7 m (0) J4/3/6 - /I 



x exp 



-/2/4/s + /i/g + /I/6 + / 3 (/ 4 2 - 4/i/e 
4/3/6 - /I 



(16) 



where 



/a = sinh- 2 [^]( 
1 



ftg 2 (0) 
2M0) 



1 



sinh 2 ^] + -A 



J\2 



2 

2 _/2>, 



+ -imu cosh^t + 0] sinh[^i](g — q' )) 



/ 2 = sinh- 2 [^]( 



2a qq (0) 

— imu cosh[0] sinh[-^t](g cosh [At] + g'sinhfAt])) 
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1 h 

fz = sinh -2 [ipt] - {imuj coshf-^t — cp] sinh^i] — sinh 2 [?/>£] — A t ) 

2 2cr gg (0) 



. , 2r , ,,-2ia gq (0)p(0) +hq(0) . l2r ,. 1 , 
/ 4 = sinh -2 [^]( 9g V Vn\ — smh 2 [^] + -£? t (g - </) 



2a qq (0) 

+ iraoj cosh[0] sinh[^t](gsinh[At] + q' cosh [At])) 
/ 5 = sinh" 2 [ipt] A t 

fe = — sinh -2 '[ipt]-(imuj cosh[ipt — 0] sinh^t] + 



2a qq (0) 



suuV 



+ A t ) 



(17) 



Here, q(t) and p(t) are the mean values of q and p, respectively, and cr qq (t), (J pp (t) and <J pq (t) 
the corresponding variances P,|T2"|| . Explicit expressions for these mean values and variances 



are given in Ref. [12]. The diagonal part of the density matrix (15) yields a Gaussian 
distribution at time t 



p(q,t) =< q\p(t)\q >= {2ua qq {t))- 1 l 2 exp[- 



1 



2a qq (t) 



{q - mf 



(18) 



where 



-Xt 



q(0) 
1 



cosh(^t) H — — — - smhiipt) 

V 



+ 



mip 



p(0) smh(ijjt) 



2m 2 X(uj 2 — Xp\ q 



+ e 



+ 



-2 At 



2d 



m (u — 2X p X)D qq — D pp — 2mX p D pq 
1 



2mu 

with the following notations 
moj 2 (X q — X p ) 



m(X q — X p ) 2mw 2 
2 l(X q -X p )C 3 + 2C^} sinh(2^) 



[(A, - X P )C 2 + 2C 3 ^] cosh(2^t) 



(19) 



Ci 



Co. 



4^ 2 

\d„ + - 



ill 



m 2 uj 2 

(Ag - Xp) 

2 uj 2 X pp moj 2 X 



1 MQH^-^Mo) 



-D 



D 



pi 



4^ 2 



A " A -(^ PP (0) - mV<7 w (0)) + WV M (0) 



+ 



^ 2 - ApA g 



2oo 2 - X n X 



ni 



^ q [D P p + m 2 u 2 D qq } 



^D P p + X 2 p muj 2 D ga + 4Xu 2 D, 



m 



'qq 



I "I 



c. 



2mip 



m uj a qq (0) + a P p(0) + 



UJ 2 - X p X q 



(X q Dpp + 2mu) D pq + m uj X p D. 



qq) 



For the values X p = X q = 0, D pp = D qq = D pq = 0, o~ pp (0) = h 2 / (Aa qq (0)) and cr qp (0) = 0, we 
obtain the same result with these expressions as in Refs. [[LJ|^(]]. For X q = 0, D pp = D qq = 
D pq = and cr pp (Q) = h 2 / (Aa qq (0)) and cr qp (Q) = 0, our results coincide in the underdamped 
limit with the results of Ref. JH]], where tunneling was studied with the inverted Caldirola- 
Kanai Hamiltonian. For X q = and D qq = 0, our result can be transformed to the one of 
Ref. |J22|. 
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III. CALCULATED RESULTS 



An influence of the friction and diffusion coefficients on the tunneling was considered in 
Refs. [|IT||12p3| . Here, we study the time-dependence of the decohering amplitude and the 



non-diagonal part of the density matrix for different sets of the transport coefficients. In 
order to demonstrate the effect of the diffusion and friction coefficients of the coordinate, 
D qq and X q , on the change of the distance between phase space trajectories, we take a simple 
expression for the diffusion coefficients: 

D pp = (1 + K)\mhuj e ff coth(hu> e ff /(2kT)) /2, 

D qq = (1 — n)h\coth(hu) e ff /(2kT))/(2mu) e ff), D pq — 0. (20) 

Here, k is an adjustable parameter. The parameter o; e // could be found from a microscopic 
consideration of the open system. With k — 1 and D qq = D pq = we obtain the "classic" 
set of diffusion coefficients which does not preserve the non- negativity of the density matrix 
at all times §,0,||. 

As an example we consider the relative motion of the two nuclei 76 Ge and 170 Er at 
the Coulomb barrier which is approximated by the inverted oscillator. Fig. 1 shows the 
time-dependence of the density matrix p(q,q') for k = and 1 in (p0|), X p = (1 + k)X 
and X q = (1 — k)A. Since with k = and 1 the density matrix is practically diagonal 
after a short time interval of about 5 x 10 -22 s, semiclassical methods work quite well in 
heavy ion collisions. The density matrix becomes faster diagonal in the case k = 1 than 
for k = 0. The time behaviour of the non-diagonal components of the density matrix 
is evidently correlated with the time-dependence of the decoherence D which is shown in 
Fig. 2. After a decrease of D during a short time interval the decoherence increases indicating 
a depression of the interference between different states (trajectories). The decoherence 
increases slowest for k = and more rapidly for higher temperatures. Further, Fig. 2 shows 
that the decoherence amplitude decreases with increasing X q . This has the consequence that 
the penetrability through the barrier increases due to a larger interference between different 
states (trajectories) |TT| , [T2"| . 

In order to show the role of D qq in distorting the coherence between states, we compare 



the time-dependence of D in Fig. 3 for k = 0, 0.5 and 1 in ( PUD with X p = 2X and X q = 0. 
For times t > 5 x 10~ 22 s, which are of interest for physical observables, the decoherence 
increases fastly for D qq — (k — 1). With D qq ^ (k < 1) the interference between different 
states survives a longer time. 

Eq. (0) can also be solved by rewriting it in a system of equations for the matrix elements 
of p in some basis ||23|| . These equations can be numerically treated for arbitrary potentials. 
With complete orthogonal set of basis functions \n > we obtain from Eq.(^) the system of 
equations for the matrix elements of density matrix p: 

= X]{-^(< rn\Ho\l > pi n - < l\H \n > p m i) 

+ PmlBln + PlnCml + Pll'A m u> n }, (21) 

I' 

where the coefficients are defined as follows 
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B ln = £(D7A+ A+ n + D?AvAj n - D 2 A w A+ n - £>+*A+ AjTj, 
v 

C ml = HflfA^ + Di*A n A ml , - £>+Aj7,A+, , - D~*A+ A",), 
/' 

An»'n = - W + ^i + )A+ z A+„ - (L>r* + Df*)A^ n + (£>j + Dj*)A^Aj7„ 

+ (Dl + Df )A ml Aj n . (22) 



Here, A~ n =< m|a|ra > and A+ n =< m|a + |n > are the matrix elements of the creation a + 
and annihilation a operators, and D x = O.B(Di ±0.5(A P — A ? )), D 2 = O.B(Di ±0.5(A P + A g )). 
For the basis related to the eigenfunctions of harmonic oscillator with the frequency u, 
D\ = (mujDqq — Dp P /muj + 2iD pq )/h and -D2 = {mujD qq + D pp /mu)/h. For calculations, 
either the eigenfunctions of harmonic oscillator or the eigenfunctions of potential U(q) are 
convenient as complete orthogonal set of basis functions \n > With the initial state of the 
open system determined by the wave function ^(q) the initial density matrix is calculated 
as p mn {t = 0) =< >< ty\n >. With this initial condition we can solve Eq . (pip and find 
the time dependence of average value F = Tr (p(t)F) of any operator F and of diagonal and 
nondiagonal elements of the density matrix. 

Let us consider a system with mass m=53m (m is the mass of a nucleon) in a symmetric 
double-well potential 

8AU 2 1QAU 4 
U{q) = -^Q 2 + -J^l 4 (23) 

with AU=1.5 MeV, L=3 fm and start with an initial Gaussian state for the density matrix 
with a variance <J qq (0) = 0.14 fm 2 in the left well at q(0) = —1.5 fm. The calculated time- 
dependence of p(q, q') is presented in Fig. 4 for k = 1 and k = in (p0|) with A p = 2A and 
X q = 0. The transition of the system to the right well mainly occurs along the direction 
q = q'. At the same time the non-diagonal part of the density matrix is larger with D qq ^ 
than with D qq = 0. For D qq 7^ 0, the distribution in the right well is wider and the transition 



rate between the two wells is larger [23,24 . 



IV. SUMMARY 

Using the path integral method for master equations of general Lindblad form for Marko- 
vian open quantum systems, we obtained an analytical expression for the propagator of the 
density matrix of a general quadratic Hamiltonian coupled linearly (in coordinate and mo- 
mentum) with the environment. The time-dependent diagonal and nondiagonal elements 
of the density matrix in coordinate representation were calculated as a function of different 
sets of transport coefficients for the inverted oscillator and double-well potentials. At times 
of interest for heavy ion collisions at the Coulomb barrier, the density matrix is practically 
diagonal which justifies the use of semiclassical methods. The time behaviour of the decoher- 
ence crucially depends on the choice of the friction and diffusion coefficients. With diffusion 
coefficients preserving the non-negativity of the density matrix at any time, the decoherence 
increases slower than in the classical case with D qq = 0. Therefore, the penetrability of a 
barrier is larger in the case of D qq 7^ due to a stronger coherence between different states. 

G.G.A. is grateful to the Alexander von Humboldt-Stiftung (Bonn) for support. This 
work was supported in part by DFG and RFBR. 



S 



REFERENCES 



[1] Lindblad G. 1976 Commun. Math. Phys. 48 119; 1976 Rep. on Math. Phys. 10 393 
[2] Belavin A. A., Zel'dovich B.Ya., Perelomov A.M. and Popov B.S. 1969 JTEP 56 264 
[3] Davies E.B. 1976 Quantum theory of open systems (New York: Academic Press) 
[4] Dodonov V.V. and Man'ko V.I. 1983 in Group Theoretical Methods in Physics, Vol. 2, 

ed. M.A.Markov (Nauka, Moscow) 
[5] Dodonov V.V. and Man'ko V.I. 1989 Reports of Physical Institute 167 (1986) 7; 191 

171 

[6] Dekker H. 1981 Phys. Rep. 80 1 

[7] Sandulescu A. and Scutaru H. 1987 Ann. Phys. (N.Y.) 173 277 

[8] Isar A., Sandulescu A., Scutaru H., Stefanescu E. and Scheid W. 1994 Intern. J. Mod. 

Phys. A 3 635; Isar A., Sandulescu A. and Scheid W. 1993 J. Math. Phys. 34 3887 
[9] Antonenko N.V., Ivanova S.P, Jolos R.V. and Scheid W. 1994 J. Phys. G: 

Nucl. Part. Phys. 20 (1994) 1447 
[10] Strunz W.T. 1997 J. Phys. A: Math. Gen. 30 4053 

[11] Adamian G.G., Antonenko N.V. and Scheid W. 1998 Phys. Lett. A 244 482 
[12] Adamian G.G., Antonenko N.V. and Scheid W. 1999 Nucl. Phys. A 645 376 
[13] Fujikawa K., Iso S., Sasaki M. and Suzuki H. 1992 Phys. Rev. Lett. 68 1093 
[14] Weiss U. 1992 Quantum dissipative systems (Singapore: World Scientific) 
[15] Razavy M. and Pimpale A. 1988 Phys. Rep. 168 305 

[16] Caldeira A.O. and Leggett A.J. 1981 1981 Phys. Rev. Lett. 46 211; 1983 Ann. Phys. 

149 374; Leggett A.J. 1984 Phys. Rev. B 30 1208 
[17] Bruinsma R. and Per Bak 1986 Phys. Rev. Lett. 56 420 
[18] Bulgac A., Dang Do G. and Kusnezov D. 1998 Phys. Rev. E 58 196 
[19] Papadopoulos G.J. 1990 J. Phys. A: Math. Gen. 23 935 
[20] Dodonov V.V. and Nikonov D.E. 1991 J. of Soviet Laser Research 12 461 
[21] Baskoutas S. and Jannussis A. 1992 J. Phys. A: Math. Gen. 25 L1299 
[22] Hofmann H. 1997 Phys. Rep. 284 137 

[23] Adamian G.G., Antonenko N.V. and Scheid W. 1999 Phys. Lett. A 260 39 
[24] Harris E.G. 1990 Phys. Rev. A 42 3685 



9 



FIGURES 



t = 0.0 




FIG. 1. Calculated time-dependence of the module of the density matrix \p(q, q') | in the inverted 
oscillator potential, reproducing the Coulomb barrier in a 76 Ge+ 170 Er collision, for k = 1 (left side) 
and (right side) in (p0|), X p = (1 + k)X and X q = (1 — k)A. The parameters are q(0) = 0, p(0) = 0, 
^a; = ?ioj e ff = 2.0 MeV, o" gg (0) = 0.7 fm 2 , m = 53mo (mo is the mass of nucleon), fi\ = 2 MeV 
and T = MeV. The initial density matrix is presented on the top. The time is given in units of 
6.582 x 10- 22 s. 
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T = MeV 




0.8 



T = 5 MeV 




0.8 

FIG. 2. Dependence of the inverse decoherence D~ l (units hT 1 ! 2 ) on time and k, used in the 
definitions of D pp and D qq in (pOj), X p = (1 + k)X and X q = (1 — k)A at T = and 5 MeV. The 
initial values are q(0) — q'(0) = 0.01 fm and p(0) = p'(0) = 0. The other parameters are the same 
as in Fig.l. The time is given in units of 6.582 x 10 -22 s. 
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FIG. 3. Time-dependence of the decoherence D (units fi 1 / 2 ) for A p = 2A, \ q = 0, T = and 
5 MeV, and k = (solid lines), 0.5 (dashed lines) and 1 (dotted lines). The initial values are 
q(0) — q'(0) = 0.01 fm and p(0) = p'(0) = 0. The other parameters are the same as in Fig.l. The 
time is given in units of 6.582 x 10 -22 s. 



12 



t = s 





FIG. 4. Calculated time-dependence of the module of the density matrix \p(q,q')\ in a dou- 
ble-well potential (see text) for k = 1 (left side) and (right side) in ( pop at T = MeV and 
X p = 2A, X q = 0. The initial Gaussian distribution (top part) with a qq (0) = 0.14 fm 2 , p(0) = 
starts with q = —1.5 fm in the left well. The further parameters are hoj e ff = 2.0 MeV, HX = 2 
MeV and m = 53mo- The time is indicated. 
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